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SUMMARY 


A  broad  analytical  study  of  the  effect  of  edge  restraint  on  the  buckling 
stress  of  geometrically  perfect,  axially  compressed  circular  cylindrical 
shells,  based  on  linear  elastic  theory,  is  presented.  An  overall  evalu¬ 
ation  of  this  effect  is  made  by  studying  the  16  limiting  cases  or  restraint 
which  correspond  to  the  vanishing  of  generalized  forces  and  generalized 
displacements  at  an  edge.  Of  the  16  limiting  cases,  7  of  which  are  studied 
for  the  first  time,  it  is  found  that  8  permit  buckling  to  occur  at  stresses 
less  then  the  classical  theoretical  value  of  buckling  stress.  The  fact 
that  both  theory  and  experiments  indicate  that  the  effect  of  shell  length 
on  the  buckling  stress  can  usually  be  disregarded  is  used  in  the  analysis. 
This  'act  suggests  that  shells  may  be  assumed  to  be  semi -infinite;  as  a 
consequence,  formulas  for  buckling  stresses  are  determined  in  essentially 
closed  form. 
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INTRODUCTION 


Recent  attention  has  been  devoted  to  the  study  of  the  effect  of  edge 
restraint  on  the  buckling  behavior  of  thin  circular  cylindrical  shells 
subjected  to  axial  compression  [1  through  6]*.  A  recent  survey  of  the 
results  of  these  studies  and  the  conclusions  drawn  from  them  has  been 
prepared  by  Hoff  [ 7 ] • 

In  addition  to  the  strict  enforcement  of  edge  conditions  at  the  instant  of 
buckling,  Stein  [8,9],  Fischer  [10],  and  Almroth  [11]  have  accounted  for 
the  prebuckling  deformations  due  to  edge  restraint  from  the  inception  of 
compressive  loading  for  shells  with  the  various  types  of  simply  supported 
and  completely  clamped  edges.  The  effect  upon  the  comp  ssive  buckling 
stress  of  lateral  pressure  was  also  studied  by  these  autnors. 

The  present  study  is  a  broad  treatment  of  the  effect  of  edge  restrairt  on 
the  compressive  buckling  stress  of  circular  cylindrical  shells;  the  analysis 
is  simplified  by  treating  semi-infinite  shells,  since  previous  experience 
has  indicated  that  the  effect  of  shell  length  on  the  buckling  stress  is 
negligible,  unless  the  shells  are  very  short  [%6].  Limiting  cases  of  edge 
restraint  which  correspond  to  the  vanishing  of  generalized  forces  and 
generalized  displacements  at  an  edge  are  treated  with  the  aid  of  the  linear 
Donnell  equations.  The  buckling  solutions  are  obtained  without  great  diffi¬ 
culty,  and  most  of  the  results  are  presented  in  essentially  closed  form. 

Of  the  16  limiting  cases  of  edge  restraint,  7  are  examined  here  for  the 
first  time;  additional  information  regarding  the  remaining  cases  is  also 
presented.  In  particular  ,  a  closed-form  expression  for  the  buckling 
criterion  treated  numerically  in  Reference  2  (and  again  in  Reference  6)  is 
derived,  the  solution  found  by  Hoff  in  Reference  1  is  extracted  from  this 
expression  as  a  special  case,  and  the  results  of  Reference  6  for  long  shells 
with  completely  clamped  edges  are  reestablished  readily. 

With  these  results  at  hand,  it  is  possible  to  summarize  the  information 
obtainable  from  linear  elastic  theory  regarding  the  effect  of  edge  restraint 
on  the  buckling  stress  of  axially  compressed  circular  cylindrical  shells. 

The  governing  equations  which  are  utilized  in  the  analysis  are  discussed 
thoroughly  in  order  to  put  the  results  presented  here  and  those  previously 
obtained  in  References  2,  8,  9/  10  >  and  11  in  proper  perspective. 

The  present  investigation  represents  a  single  detailed  analysis  ir  the 
overall  program  of  plate  and  shell  theoretical  analyses  currently  in 
progress  at  Stamford  U.  versity  in  the  Department  of  Aeronautics  and  Astro¬ 
nautics.  Unfortunately,  there  appears  to  be  no  published  experimental  data 
available  at  the  present  time  with  which  the  analytical  results  presented 
herein  can  be  compared.  A  systematic,  complementary,  experimental  investi¬ 
gation  into  edge-fixity  effects  is  currently  in  progress  under  this  same 


Numbers  in  brackets  indicate  the  references  listed  at  the  end  of  the 
report . 
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contract,  but  quantitative  results  are,  as  yet,  unavailable.  In  view  of 
the  success  of  the  statistical  approach  to  shell  buckling  experiments 
developed  earlier  in  the  experimental  phases  of  the  contract  work,  reason¬ 
able  correlation  of  test  results  with  linear  elastic  theory  can  be  expected. 


BASIC  THEORY 

THE  KARMAN-DONNELL  EQUATIONS 


The  plate  equations  of  von  K^rm^n  were  extended  by  Donnell  [12]  to  the 
analysis  of  circular  cylindrical  shells;  these  equations  are  the  foundation 
of  the  present  analysis  and  may  be  written  as 


and 


E 
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W  ' 
,xx  ,yy 


(1) 


(r)v4W  -  F  W  +2F  W  -F  (W  +i)=0 
t'  ,yy  , xx  ,xy  ,xy  ,xx  ,yy  R 


where  F  is  an  Airy  stress  function  defined  such  that 


(2) 


,xx 


(3) 


Boundary  conditions  for  complete  cylinders  with  edges  along  constant  x- 
curves  are  defined  by  the  four  following  alternate  choices  of  quantities 
which  may  be  prescribed: 
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0 

T 

*y 

R  =  -  D(W  +  (2-v)W  )  +  to-?W  +  tT°  W 

x  ,xxx  ,xyy  x  ,x  xy  ,y 
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or 


M  =  -  D(W  +  vW  ) 
x  ,xx  ,yy 


With  reference  to  the  figure  (page  3),  U,  V,  and  W  are  the  components  of 
the  displacement  vector  in  the  axial,  tangential,  and  radial  directions, 
respectively,  and  the  coordinates  describing  the  undeformed  shell  midsur- 
ffice  are  x,  y  =  R 9,  and  z. 


The  strains  at  the  shell  midsurface  are  given  by 

4  --  U,X  +  5<w,/ 


(continued) 
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=  V  -|  +  i(W)2( 
,y  R  2  '  ,y  ) 
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The  changes  in  curvature  and  twist  of  the  shell  midsurface  are 


k.  =  -  W 


<  =  -  W 


X  ,xx  '  y  ,yy 


**y  ■  -  “,xy 


The  stress-strain  and  moment -curvature  relations  for  an  isotropic  Hooke  an 
elastic  material  are 
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where  Mx  ,  are  the  moment  resultants  and  D  =  EtJ/l2(l-v  )  . 


THE  LINEARIZED  DONNELL  EQUATIONS  FOR  COMiftESSIVE  AXIAL  PRESTRESS 


The  membrane  solution  to  the  K  arman  -  Donnell  equations  for  uniform  axial 
compression  is 

*<,  --4 


w  _  vc£ 

0  E 

U  =  -  2* 

0  E 


v0  =  0 


(8) 


This  membrane  or  "momentless"  state  can  exist  only  if  the  uniform  radial 
expansion  WQ  is  permitted  to  occur. 


4 


1 


If  equilibrium  states  vhich  differ  only  slightly  from  the  membrane  state 
are  considered  and  if  sufficiently  small  changes  in  all  quantities#  say# 

F,  #  U,#  Vi#  and  W,#  are  assumed#  the  basic  equations  can  be  linearized, 
lfie  substitution  or  Fq  +  F,#  etc.#  into  equation  (l)  through  (7)  and 
subsequent  linearization  result  in  the  following  equations  governing  the 
perturbations 'from  the  membrane  state: 


2  w 

R  "1#  xx 


(9) 


and 


+  ow. 


l#xx 


A 


xx 


(10) 


where  the  increments  in  the  midsurface  stresses  are  obtained  from  F^  by 
the  relations 


(11) 


These  stresses  are  related  to  the  displacement  components  by  the  following 
equations : 


dWui,y  +  vi#x)  "  J 

Equations  (9)  through  (12)  have  been  put  in  a  convenient  nondimenalonal 
form  by  Nachbar  and  Hoff  [2]#  who  define  the  following  new  variables*: 


Hie  notation  used  here  is  slightly  different  from  the  notation  in 
reference  2. 
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where  orc^  Is  the  classical  buckling  stress  from  linear  theory  given  by 
the  formula 


cl 


[3(1-/)] 


_  t 

172  R 


(in) 


With  the  aid  of  the  definitions  in  equation  (13),  equations  (9)  through 
(12)  can  be  replaced  by  the  following  dimensionless  expressions: 


A  =0 

(15) 

V^w  +  2pw  -  f  ..  =0 

>55  >55 

(16) 

0  0 

Ef  >  T  ■=  -Ef  t  ,  a  =  Ef 

,u 

(17) 

>99  xy±  >  59  yx 

u  -  «  f  _•  vf  „  ^ 
>5  >99  >55 

>9  >55  >99 

u  m  +  v  -  -2(l+v)f 
>9  >5  >59  J 

> 

(18) 

where  p2  (  )  *  (  )  ««  +  (  )  __ 

>55  >W 


(19) 


Hie  boundary  conditions  for  the  linearized  equations  require  the  specifi¬ 
cation  of 


and 


u  or  a  ■  Ef 
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=  -Et(Si)  [*,ill+(2-v)\lW+2p\i) 


M 

x 


+  vw 


>99 


(20) 


Equations  (15 )  through  (18),  (19)>  and  (20)  will  be  used  in  the  subsequent 
buckling  analysis;  they  are  equivalent  to  the  equations  deduced  by  Donnell 
[13]  for  compressive  axial  presiress. 

The  purpose  of  describing  in  detail  the  derivation  of  the  basic  equations 
is  to  emphasize  the  following  two  points: 
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1.  Consistent  linear  stability  equations,  along  with  their  assoc¬ 
iated  boundary  conditions,  are  always  obtained  by  linearizing 
consistent  nonlinear  equations  about  the  basic  state  of  equili¬ 
brium  whose  ^ability  is  to  be  investigated.  The  alternate 
techniques  for  deriving  stability  equations,  even  when  based 
upon  appropriate  modification  of  small  displacement  theory  by 
the  inclusion  of  the  effects  of  "apparent  buckling  forces",  must 
be  utilised  with  caution  and  are  not  recommended;  the  accuracy 
and  consistency  of  the  resulting  equations  and  associated 
boundary  conditions  cannot  be  assured  or  easily  assessed.  The 
boundary  conditions  employed  in  Reference  2,  for  example,  are 
inconsistent  with  the  differential  equations  upon  which  the 
analysis  is  based,*  as  has  been  pointed  out  in  Reference  6;  a 
closed-form  expression  for  the  buckling  criterion  studied  in 
Reference  2  will  be  derived  subsequently,  thus  demonstrating  the 
advantages  of  consistency. 

2.  The  linear  stability  equations  used  by  Stein  [8,9],  Fischer  [10], 
and  Almroth  [11]  in  their  studies  of  compressive  buckling  of 
cylindrical  shells  were  derived  by  linearization  of  the  Kirmin- 
Donnell  equations  about  basic  states  of  equilibrium  which  differ 
from  the  membrane  state  equation  (8).  These  authors  account  for  the 
effect  of  axl symmetric  prebuckling  deformations,  as  previously  men¬ 
tioned  in  the  introduction,  Which  rigorously  enforce  the  prescribed 
boundary  conditions  in  the  basic  states  of  equilibrium.  For 
axially  compressed  cylinders,  Almroth  [11],  who  also  performed 
calculations  with  these  nonuniform  prebuckling  deformations 
neglected,  states  that  the  results  which  he  obtained  using  the 
equivalent  of  equations  (9)  and  (10)  are  not  significantly 
different  from  those  obtained  using  equations  in  which  the  bound¬ 
ary  conditions  are  rigorously  enforced  from  the  inception  of 
compressive  loading. 


SOLUTIONS  OF  THE  BASIC  EQUATIONS  FCB  COMPLETE  CYLINDRICAL  SHELLS 
THE  FORM  OF  THE  SOLUTIONS 

In  the  case  of  the  complete  cylindrical  shell,  f  and  w  must  be  periodic 
functions  of  the  circumferential  coordinate  y.  Since  only  derivatives  of 
even  orders  appear  in  equations  (15)  and  (16),  it  is  obvious  that  products 
of  functions  of  {  alone  with  either  sin(sy/R)  or  cos(sy/R),  where  s 
is  an  Integer,  can  be  made  to  satisfy  them.  The  ordinary  differential 


Although  the  problem  was  formulated  with  inconsistent  boundary  conditions, 
the  numerical  calculations  were  performed  after  the  inconsistent  terms 
had  been  discarded  on  the  basis  of  an  order-of -magnitude  argument. 
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equal  ions  for  functions  of  I  wi  Jch  ore  ottoia*,ri  l >y  aurh  an  08C.oed 
colution  are  solved  I  y  exponential  nmetlonu  cl  \i  thus,  without  loss  or 
generolity,  solutions  are  cujbuiim  d  of  the  1'orm 


p* 

w  t  Ac  cob  i# 

pi 

1’  *  Ikr  cob  np 


(m) 


where  A  and  B  are  constants;  and,  free  equation  (13),  n  oust  Ur 
related  to  a  by  the  equation 

n  “  B{'2r)  ‘  •(12(l-v2))‘l/‘*(i)1/c  (22) 

By  introducing  equations  (21^  into  equations  (15)  and  (16),  the  ho*K>geneous 
equations  for  the  values  of  p  are  found  to  be 


-  2“ 

P2  (p2-n2) 

-  - 

A 

((p2-n2)2  ♦  2pp2]  -p2 

—  • 

B 

*  0 


(23) 


For  nontrivial  solutions  to  equations  (23) .  the  deterainant  cf  the  matrix 
of  coefficients  suet  vanish,  ifclch  gives  the  characteristic  equation 


(p2-r2)  I(p2-n2)  ♦  2pp2  J  ♦  p*  •  0  (2-) 

2  2 

If  p  =  n  ,  then  only  the  trivial  solution  results;  by  excludlr^;  this 
possibility,  equation  (2>t)  con  be  written  as 

C2  ♦  2pC  ♦  1  .  0  U‘j) 


where 


,  2  2\ 

c.l&p- 

p 


From  equations  (23)  it  is  found  that 


(26) 


(27) 


where  C  Is  a  solution  to  equation  (25). 
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THE  ROOTS  or  RUATIOM  (24)  for  p  <  1 

Solutions  to  the  characteristic  equation  (2k)  were  first  given  by  Nachbar 
[H*J  without  the  details  of  their  derivation.  A  concise  discussion  of  the 
solutions  will  be  given  here*  for  the  purpose  of  providing  various 
Identities  that  will  be  used  subsequently. 


If  pal, 


the  two  solutions  of  equation  (25) 

p  1/2 

cx  -  -  P  ♦  i(l-P  ) 

2  1/2 

c2  -  P  -  i(l-P*) 


can  be  conveniently  written 


s  Q 

(28) 

=  Q 

(29) 

where  Q  and  Q  are  complex  conjugates  and  i  -  +  \T- T.  The  modulus  of 
Q  is  unity;  that  is, 


Oft  =  1  (30) 

The  following  expressions  for  the  positive  square  roots  of  Q  and  Q  are 
noted: 


Q 1/2  =  +  -^-[(1-P)1/2  +  i(l+p)l/2] 
42 


(31) 


Ql/2  -  +-r-[(l-p)l/2  -  i(l+p)l/2 J  (32) 

42 


FSpob  aquations  (26),  (26)  and  (29),  four  quadratic  equations  for  the  roots 

p  are  obtained;  that  is 


2  2 
E.  ■-»-  = 
P 


:  «1/2 

♦ 

.  «1/2 

+ 


(33) 

(34) 

(35) 

(36) 


The  two  roots  of  equation  (13),  designated  p^  and  ,  are 

Pl  -  |  [-«1/2  -t^n2)l/2]  (37) 

The  notation  employed  here  differs  from  that  of  Reference  14;  the  notation 
of  Reference  4  is  utilized. 


9 


(38) 


P5  -  |  [-«1/2  *  (Q^n2)^] 

Similarly,  the  remaining  six  roots  are  obtained  from  equations  (34),  (35), 
and  (36).  They  are 


from  equation  (34) : 

~  P2 

=  -P5 

-  p6 

II 

1 

from  equation  (35): 

P3 

P1 

L  P? 

=  p5 

from  equation  (36): 

r  14, 

=  -  P5 

L_  Pf, 

”  -  pl 

2  1/a 

The  quantity  (Q+4ir)  is  taken 

(Q+4n2)  and  may  be  written  in  the 

to  be  the  positive  square  root  of 
following  way: 

(Q+4n2)  =  -^-[(r+4n2  -p)^2  +  i(r-4n2  +  p (40) 

4  2 

where 

r  s  +(l-8n2p  +  16^)  =  +[(4n2-p2)  +  1-p2]  (4l) 

Note  that  if  p  s  1,  then  r  must  be  positive. 

With  the  aid  of  equations  (31)  and  (40),  p.  and  p2  can  be  written  with 
their  real  and  imaginary  parts  explicitly  separated;  the  expressions  are 


and 


where 


-Mai +  iei’ 

2  V2 

(42) 

-V(02  +  1P2> 

2  V2 

(43) 

“l  = 


.(l-p)1/2  -(r+4n2-p )1//2;  ^  =  -(l+p)1//2  -(r-4n2+p)1//2 
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a2  .  (l-p)1/2  -(r+n2-p)1//2]  &2  =  (l+p)1^  -(r-brf+p)1^  (44) 

The  following  identity  is  useful  and  may  be  verified  by  direct  expansion: 

(r+4n2-p)1//2(r-4n2+p)1//2  =  (1-p2)1^2  (45) 

SOLUTIONS  SUITABLE  FOR  ANALYZING  THE  BEHAVIOR  OF  LONG  SHELLS 


Four  roots  of  equation  (24),  p  through  p,  ,  have  negative  real  parts 
for  values  of  p  less  than  unity;  a  proof  or  this  statement  is  given  in 
Appendix  I.  Roots  of  this  type  indicate  that  solutions  which  decay 
exponentially  with  f  exist  and  that  sufficiently  long  shells  may  be 
treated  mathematically  as  semi -infinite  whenever  p  <  1.  A  discussion  of 
the  applicability  of  "semi -infinite”  buckling  results  to  shells  of  finite 
length  is  given  in  the  next  section. 


If  a  shell  is  assumed  to  be  semi -infinite  along  the  positive  x-axis  with 
its  edge  at  x  -  0,  f  and  w  will  be  bounded  as  x  —  ®  and  can  be 
written,  with  the  aid  of  equations  (27)  through  (30),  (37)#  and  (35),  as 
follows  when  psl: 


where 

into 


w  -  l^e^l*  +  +  AjC^l^  +  A^e^2^j  cos  ncp 

f  -  -  j^Q'Vl5  +  AgCf^V  +  AjQ  ‘Vl1  +  A4Q  ■V2IJ 

=  -  Ja^I6  +  A +  AjQe^l^  +  A4Qep2*J  cos  nqp  (47) 

A.  through  A  are  arbitrary  constants.  Substituting  equation  (47) 
eqifttion  (17)  yifelds  the  following  results: 


(46) 


cos  np 


cos  ncp 


1  ^5ePl5  + 

Ag5eP2S  +  AjQePl5  +  A4QeP2^  J  cos  ncp 

(48) 

n 

;  +  A2p2QeP25  +  A3P1QePl6  +  A4P2QeP2^ 

sin  nqp 

L 

\^9) 

r"iLp^QePlt  i 

i  A2P^QeP2e  +  A^QeV-  +  A^QeVl 

cos  rKp 

L 

J 

(lv0) 

II 


Substituting  equations  (46),  (48) .  and  (‘jO)  into  the  expression  for  v 
in  equations  (l8)  gives 


v  -  f  . .  -  vf  +  w 
>«P  ,H 


A1[l-Q(p^+vr^)]ePl^  +  A?  [l-Q^+vn1)  ]ep2 

+  A3[l-Q(p^  +  vn‘  )]ePl^  +  A^  [1-Q(p2  +  vn2)]J  cos  ncp  (5l) 

The  utilization  of  equations  (48)  and  (50)  in  the  expression  for  u  ,  , 
equations  (18),  and  subsequent  integration  of  the  result  yields 


7T/  2A,  2x,  p06 


u  = 


-  p  F  -  (n2+VP^  p  | 

A  o  - ±-  eV  +  A  Q  - —  ep25 

1  Pt  2  pr 


■1 

2.  -2 


2 

2.. -2, 


(n  +Vpi}  p  |  (^+^2)  ~  ? 

+  A  Q  - —  ePl£  +  A,  Q,  - —  eP25 


Po 


cos  ncp 


(52) 


The  process  of  integration  with  respect  to  £  could  introduce  a  linear 
function  of  |  in  addition  to  the  terms  given  in  equation  (52).  This 
linear  function,  however,  represents  the  same  form  of  displacement  as  Uq, 
and  thus  it  is  uniquely  included  in  Uq  when  the  axial  stress  a  is 
fixed  in  magnitude. 

Equation  (52)  can  be  put  in  a  more  convenient  form  by  the  use  of  the 
identity 


P1P2 


n 


(53) 


which  may  be  verified  by  direct  expansion.  By  using  this  result,  equation 
(52)  can  be  written  as 

u  =  [  A1Q(p2+vp1)ePl^  +  A2Q(p1+vp2)eP2l 

+  A3Q(p2+vp1).  Plfj  +  A^QXp^+vp^e1^  ]  cos  ncp  ( ) 

All  the  necessary  physical  quantities  can  new  be  obtained  from  equation:: 
(46),  (48)  through  (5l);  and  (54)  for  the  ueml-infirii  t.e  shell  in  term n  of 

the  four  arbitrary  constants  A.(  through  . 
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APPLICABILITY  OF  SEMI-INFINITE  SHELL  BUCKLING  SOLUTIONS  TO  SHELLS  OF  FINITE 
LENGTH  - - 


In  the  preceding  discussion  of  solutions  to  the  governing  differential 
equations,  the  solution  suitable  for  the  analysis  of  semi -inf inite  shells 
was  assumed  to  he  adequate  for  "long"  shells.  A  shell  is  "long"  according 
to  this  definition  if  the  conditions  at  one  of  its  edges  can  be  analyzed 
without  regard  to  the  conditions  prevailiig  at  the  other. 

It  is  well  known  that  there  are  two  decay  lenghts  for  edge  disturbances  in 
circular  cylindrical  shells;  these  correspond  to  the  exponentials  propor¬ 
tional  to  Oj_  and  CK^.  Since  the  quantity  n  is  substantially  less  than 
unity  for  thin  shells,  it  follows  from  the  proof  in  Appendix  I  that  the 
-  exponential  terms  in  the  semi -infinite  shell  solution  have  large  decay 
lengths,  thus  suggestiig  that  shells  may  have  to  be  very  long  before  this 
solution  alone  can  satisfactorily  describe  their  buckling  behavior. 

Results  obtained  in  References  5  and  6,  however,  indicate  that  buckling 
stresses  are  essentially  independent  of  shell  length  for  shells  of  practical 
geometric  proportions  in  the  cases  of  edge  restraint  that  were  investi¬ 
gated*;  buckling  stresses  obtained  in  References  2  and  4  for  semi-infinite 
shells  agree  well  with  their  finite-length  counterparts  for  all  but  the 
very  short  class  of  shells.  The  buckling  mode  shapes  in  some  cases  offer 
intuitive  Justification  for  this  agreement  between  "finite-length"  and 
"semi -inf inite"  buckling  stresses.  The  modes  for  finite -length  shells 
found  in  Reference  5,  for  example,  can  be  visualized  as  having  been  con¬ 
structed  by  appropriately  Joining,  in  a  continuous  manner,  two  oppositely 
directed,  semi-infinite  modes  at  the  midlength  cross  section  of  the 
cylinder;  the  superposition  of  a  rigid  body  rotation  and  a  reflection  of 
one  of  the  semi-infinite  modes  could  represent  the  finite-length  mode  which 
is  antisymmetric  with  respect  to  the  midlength  cross  section. 

Since  past  experience  points  to  the  conclusion  that  semi-infinite  shell 
analyses,  which  offer  great  mathematical  simplicity,  predict  nearly  the 
same  buckling  stresses  as  their  more  involved  finite-length  counteiparts, 
the  semi-infinite  assumption  will  be  used  exclusively  in  the  subsequent 
determination  of  buckling  stresses.  Little  generality  is  sacrificed  in 
adopting  this  assumption,  since  very  short  shells,  the  class  of  shells  for 
which  this  assunption  has  been  shown  to  be  in  error  in  Reference  6,  are 
likely  to  buckle  well  into  the  plastic  range. 


There  is  one  exceptional  case  that  was  found  in  Reference  6.  This  is  the 
case  of  a  relatively  short,  completely  unrestrained  cylindrical  shell, 
and  it  will  be  discussed  in  a  subsequent  section. 
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CLASSIFICATION  OF  EDGE  CONDITIONS 


SCOPE  OF  THE  INVESTIGATION 


When  cylindrical  shells  are  utilized  as  structural  components  of  flight 
vehicles,  their  edges  are  restrained  by  adjacent  structural  components, 
which  may  be  represented  by  "effective"  elastic  or  inelastic  springs.  When 
cylindrical  shells  are  tested  in  compression  under  laboratory  conditions, 
loads  may  also  be  supplied  to  the  shell  edges  by  friction  between  the 
source  of  load  application  and  the  shell  or  its  mounting  fixtures.  Such 
effects  will  not  be  considered  here;  the  types  of  edge  fixity  to  be  treated 
may  be  regarded  as  limiting  cases  of  more  realistic  edge-fixity  conditions. 

The  edge  conditions  characterized  by  equations  (20)  can  be  conveniently 
separated  into  two  categories  which,  for  the  lack  of  better  terminology, 
will  be  referred  to  as  "normal"  and  "tangential"  edge  conditions.  Normal 
edge  conditions  refer  to  the  conditions  involving  quantities  that  govern 
the  motion  of  the  shell  idge  in  the  direction  normal  to  the  tangent  plane 
of  the  shell  midsurface  at  any  point;  these  are  then  .conditions  describing 
the  behavior  of  an  edge  due  to  bending  and  twisting.  Tangential  edge 
conditions,  on  the  other  hand,  govern  the  motion  of  the  shell  edge  in  the 
tangent  plane  to  the  shell  midsurface  at  any  point  and  characterize  the 
behavior  of  the  edge  under  loads  parallel  to  the  tangent  plane.  Since 
only  limiting  cases  of  practical  edge  restraint  are  to  be  studied  here, 
these  two  categories  will  be  classified  separately < 

CLASSIFICATION  OF  NORMAL  EDGE  CONDITIONS 


Each  of  the  following  four  classes  of  normal  edge  conditions  will  be 
identified  by  two  letters  which  are  abbreviations  for  the  terminology 
assigned  to  the  particular  class: 

FF:  The  "free-free"  class  for  which 

R  =  0  and  M  =  0  (55) 

X  X 

FC:  tile  "free -clamped"  class  for  which 

R  =  0  and  w  ,  =  0  (56) 

x 

SS:  The  "simply  supported"  class  for  which 

w  =  0  and  Mx  =  0  (57) 
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* 


CC:  The  "completely  clamped"*  class  for  which 

w  =  0  and  w  =  0  (58) 

It  ahould  be  mentioned  that  within  the  linearized  theory  of  static  stability, 
the  condition  w  ■  0  in  equations  (57)  and  (58)  may  be  replaced  by 

W  =  WQ  +  W]_  =  0  (59) 

If  the  perturbation  displacement  W.  is  chosen  to  nullify  the  uniform 
radial  expansion  WQ,  then  the  displacement  W.  will  approach  infinity  at 
the  buckling  stress.  Each  of  these  equivalent  ways  of  prescribing  the 
boundary  conditions  leads  to  the  same  unique  buckling  stress,  since  the 
buckling  condition  is  defined  by  the  vanishing  of  the  same  determinant  in 
both  instances. 

CLASSIFICATION  OF  TANGENTIAL  EDGE  CONDITIONS 


Each  of  the  four  classes  of  tangential  edge  conditions  will  be  identified 
by  an  integer;  the  selection  of  one  class  of  normal  edge  conditions  to¬ 
gether  with  one  of  the  following  classes  of  tangential  edge  conditions 
represents  a  complete  description  of  edge  restraint: 


1: 

0 

cr 

=  0 

and 

0 

T 

=  0 

(60) 

X1 

2: 

u 

=  0 

and 

0 

T 

=  0 

(61) 

^1 

3: 

c° 

=  0 

and 

V 

=  0 

(62) 

X1 

4: 

u 

=  0 

and 

V 

=  0 

(63) 

As  an  example  e  classical  buckling  stress  is  obtained  for  SS3  edge 
restraint  when  equations  (57)  and  (62)  are  enforced  at  the  ends  of  the 
shell. 


"Clamped"  is  a  conventional  description  of  this  class;  the  term  "com¬ 
pletely  clamped"  is  introduced  to  provide  a  contrast  with  the  "free- 
clamped"  definition  given  previously.  No  standard  terminology  has  been 
agreed  upon,  however;  the  authors  of  Reference  6,  for  example,  refer  to 
this  class  of  restraint  as  "rigid-fixed". 
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BUCKLING  STRESS  RATIOS  FOR  LIMITING  CASES  OF  EDGE  RESTRAINT 
FF1  CASE 

This  case  was  treated  initially  by  Hoff  [l]  for  axisymmetric  buckling 
deformations.  The  generalization  to  nonaxisymmetric  buckling  wus  subse¬ 
quently  accomplished  by  Nachbar  and  Hoff  [2];  in  this  study,  however, 
boundary  conditions  which  were  inconsistent  with  the  differential  equations 
that  the  authors  utilized  in  the  analysis  were  employed  to  formulate  the 
problem.  This  inconsistency  was  later  pointed  out  by  Hoff  and  Soong  [6] 
and  was  corrected  in  their  extended  analysis  which  accounted  for  the  effect 
of  finite  shell  length. 

Although  the  correct  numerical  result  for  long  shells  was  essentially 
obtained  in  Reference  2  and  was  verified  in  Reference  6,  it  is  of  theoret¬ 
ical  interest  to  reinvestigate  this  case,  for  it  will  subsequently  be  shown 
that  a  closed-form  expression  for  the  buckling  criterion  is  obtainable. 
However,  the  important  result  of  the  analysis  is  the  procedure  of  simplifi¬ 
cation  which  can  be  used  directly  to  treat  other  cases  without  great  diffi¬ 
culty.  A  collection  of  identities  which  will  be  used  in  the  following 
analysis  appears  in  Appendix  II. 

With  the  aid  of  equations  (46),  (48),  (49),  and  (20),  the  requirements 
given  in  equations  (55)  ar.d  (60)  for  the  edge  at  x  (or  5)  =  0  lead  to 
the  following  four  algebraic  equations : 
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T^e  rog 


x,  '  Txy 


s  of  the  matrix  of  coefficient  correepond 
,  Mx;  and  ?x,  respectively. 


to  the  vanishing  of 


Equation  (64)  represents  a  system  of  homogeneous ,  linear  algebraic 
equations;  nontrivial  solutions  for  the  constants  through  A^  are 
possible  only  if  the  determinant  of  the  matrix  of  coefficients  vanishes. 
The  value  of  the  stress  ratio  p  for  which  the  determinant  vanishes  are 
critical  values  which  correspond  to  buckling  within  the  scope  of  linear 
static  stability  theory. 


Since  the  vanishing  of  the  determinant  of  the  mr.trix  of  coefficients  i6  the 
only  factor  of  interest  here,  it  it  desirable  to  multiply  the  rows  and 
columns  of  the  determinant  by  nonzero  quantities  in  order  to  obtain  alge¬ 
braic  simplifications.  The  first  operations  to  be  performed  are: 

1.  Multiply  the  columns  of  the  matrix  by  pg,  p^,  p^,  and  p^, 
respectively;  and  in  view  of  equation  (53),  n^  ractors  may  be  dis¬ 
regarded  in  the  last  three  rows  of  the  matrix  since  they  do  not 
influence  the  vanishing  of  the  determinant. 

2.  Subtract  the  first  column  from  the  second,  and  subtract  the  third 
column  from  the  fourth. 


The  modified  determinantal  equation  resulting  from  the  above  operations  is 


P2ft  Qfo-L-Pg)  P2Q 

§  0  Q 


Q(Pl“P2) 


0 


=  0 


(Px-vp2)  (l+v)(p2-P1)  (Pj."^)  (I+vK^-P-l) 


p2+[2p-(2-v)n2] 

,  2  _2\ 
(p2-pi) 

p2+[2p-(2-v)n2] 

/— 2  —2  \ 
P2~P1 

(65) 

The  factor  (pg-p-^) 
removed,  since  they 
by 

multiplying  column  2  and  the  factor  (p2-P1)  may  be 
represent  multiplication  of  the  remaining  determinant 

(p2-pi 

)(P2-P1)  =  1 

(66) 

Thus,  equation  (65) 

can  be  written 

as 

p* 

-Q 

p2Q 

-Q 

Q 

0 

Q 

0 

(PrVT2) 

(l'+v) 

(Pl-^2) 

(1+v) 

=  0 

p2+[2p-(2-v)n2] 

(px+p2) 

p2+[2p-(2-v)n2] 

(Pj^) 

(67) 
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The  following  operations  are  now  performed: 


1.  Multiply  column  2  by  p ^  and  add  the  result  to  column  1;  multiply 
column  4  by  p^  and  add  the  result  to  column  3. 

2.  Equation  (53)  and  the  following  identities  are  used  to  simplify  the 
modified  first  and  third  columns: 


Px  +  P2  =  2n  +  Q 


-2-2  2  - 

Pi  +  Pg  =  2n  +  Q 


The  determinants!  equation  becomes 


0 

Q 

W 

Q  +  k 


-Q 

0 

(1+v) 

(Pl+P2) 


0 

Q 

(pi^2) 
Q  +  k 


where 


k  =  2p  +  (l+v)nc 


-Q 

0 

(1+v) 


=  0 


(68) 

(69) 


(70) 


(71) 


If  the  rows  and  columns  of  the  determinant  in  equations  (70)  are  rearranged, 
the  result  can  be  conveniently  written  as  the  sum  of  two  determinants 


A  =  &i  +  Ag  =  0 


(72) 


where 


*1  " 


and 


Q  Q  0  0 

Q  Q  0  0 

0  0  -Q  -Q 

(P-L+Pg)  (p-l+Pj)  (1+v)  (1+v) 


(l+v)(Q-Q)2(ftfQ) 


(73) 
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Q 

Q 

0 

0 

0 

0 

-Q 

-Q 

*2  = 

(P-L+Pg) 

(V*2} 

(1+v) 

(1+v) 

k 

k 

(pi*52) 

(P!+P2) 

=  k(l+v)(Q-Q)2-[(i1+i2)Q-(p1+P2)Q]2 

(74) 


The  substitution  of  equations  (73)  and  (74)  into  equation  (72)  yields  the 
simple  result 

(ft-Q)2(l+v)[(Q-^)  +  k]  =  [(p-L+Pg)^  -  (p:+P2)Q32  (75) 

However, 

Q  +  Q  +  k  =  (l+v)n2  (76) 

so  equation  (75)  offers  the  two  possible  solutions 

(Q-Q)(l+v)n  =  iUp-j+p^Q  -  (P1+P2)Q.]  (77) 

By  direct  expansion,  the  following  expressions  are  easily  obtained: 

(Pl-*P2)5  -  (p1'»P2)Q  =  -  ~^[p(r-4n2+p)1^2  -(l-p2)1^2(r+4n2-p)1^2]  (78) 

(Q-Q)  =  -  2i(l-p2)l/2  (79) 

Equation  (77)  is*  therefore,  equivalent  to 

(1-p2)1/2  ^2(l+v)n  =  ±[p(r-4n2+p)1^2  -  (l-p2)1^2  (imJu^-p)1/2]  (80) 

2  l/2 

The  quantity  (r+4n  -p)  '  cannot  vanish  unless  p  =  1.  It  may  also  be 
observed  that  p  =  1  will  satisfy  equation  (80).  Since  solutions  for  p  <  1 
are  of  interest,  the  multiplication  of  equation  (80)  by  (r+kn^-p)1/2  is  per¬ 
missible;  and,  if  the  common  factor  (l-p2)l/2  is  i  nored  after  equation  (45' 
is  employed,  the  following  two  alternatives  are  deduced  from  equation  (80) : 

V2  (l+v)n(r+4n2-p)‘l,//2  =  p-(r+4n2-p)  (8l) 

*J~2  (l+v)n(r+4n2-p)1//2  =  -p+(r+4n2-p)  (82) 

The  smallest  value  of  p  is  obtained  from  equation  (82);  it  was  found 
numerically  in  Reference  6  to  be  0.380  for  v  =  0.3  and  corresponds  to 
n  =  0.30. 


19 


The  above  equations  orr  txtrmly  oenaitive  to  epproxleate  techniques  of 
solution.  However.  If  n  it  pencilled  to  vanish  In  either  equation  (81) 
or  (KP',  which  lap  Lies  axiryaaetrlc  deforaatlon,  the  solution  p  •  1/2 
results;  this  la  the  solution  that  vu  firat  obtained  by  Hoff  in  Reference  1. 

In  their  study  of  finite -length  effects,  Hoff  and  boong  [6]  discovered 
that  very  short  shells  with  FT1  edges  car.  be  buckled  by  extreaely  low 
values  of  stress;  the  aode  of  instability  Is  one  of  twisting  as  a  curved 
strip  and  can  easily  be  visualised.  This  aode  of  buckling,  however,  should 
not  be  a  Better  of  concern  for  thells  that  are  currently  of  Interest. 

FTC  CASE 

Equations  (55)  and  (61 )  will  be  satisfied  at  |  -  0  if  the  first  row  of 
the  coefficient  aatrlx  In  equations  (6k)  is  replaced  by 

IppW  Pji  p2Q  PjQ]  (83) 


If  the  saae  sequence  of  operations  Is  performed  which  slspllfled  the  Krl 
case  to  the  fors  given  lr.  equations  (70),  the  buckllrg  condition  for  an 
FT2 -supported  shell  car.  be  written 


-5n* 

-5(pi*p2) 

■Qn2 

“Q(pl  *Pg  ) 

Q 

0 

Q 

0 

(PX-P2) 

(l*v) 

(pi-V 

(l*v) 

Q  ♦  k 

5  k 

(Pj^) 

p 

If  row  2  Is  sultlplled  by  n  and  added  to  row  1  In  equation  (8b)  and  the 
rows  and  coluens  are  rearranged,  the  buckling  condition  can  be  written 


A  *  »  0 


(85) 


where 


$  k  0 

Q  5  o 

o  0  -^(p1*p2) 

(p^)  (pa*p2)  d*y) 


0 


0 

-5(Pi*p2) 

(l*v) 


-(l*v)tf-Q2)[Q(p1^P2) 
■Q(p^*p2  ) ) 

(86) 


and 


t 

] 
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or 


Q 

0 

(p-j+p^) 

k 


Q 

0 

W 

k 


0 

-Qfip-j+Pg) 
( 1+v ) 

(P1+P;J) 


0 


-Q^+p^) 


(1+v) 


(pl+p2 


•^2  =  (Q“Q)k(l+v)  [Q(p1+p2)-;dP1+P2)  ^"(^“^[^(Pj+P^-QXp^Pg)  3  (87) 

by  use  of  the  identity 

(P1+P2)(P1+P2)  =  r  (68) 

'Fie  substitution  of  equations  (86)  and  (87)  into  (85),  along  with  the  use 
of  equation  (76),  gives 

(Q-Q)  {(l+v)2n2[Q(p1+p2)-Q(pi+p2)]-r[Q(p:L+p2).Q(p1+p2)]}  =  0  (pc,) 

The  quantity  [ (Pt+P2)Q-(p-,+P2)Q3  is  given  in  equation  (78),  and  direct 
expansion  of  the  other  terms  leads  to 

[Q(P1+P2)-Q(P1+P?)]  =  [p(r-4n2+p)1/2+(l-p2)1/2(r4iin2-p)1//2]  (90) 

V2 

If  equation  (89)  is  multiplied  by  (r+4n2-p)1//2  and  constant  terms  and  fac¬ 
tors  which  vanish  at  p  =  1  are  disregarded,  the  buckling  condition  becomes 

[r+(l+v)2n2]p  =  [r-(l+v)2n2](r+4n2-p)  (91) 

or 

r2  +[(4-(l+v)2)n2  -  2p]r  -  4(l+v)2ni+  =  0  ( 92 ) 

The  quantity  n  is  proportional  to  t/R,  and^thus  it  is  small  corn}  .red 
to  unity  for  thin  shells.  By  neglecting  the  11  term  in  equation  (92)  and 
by  noting  that  r  can  vanish  only  for  p  >  1  ,  an  approximate  solution  to 
equation  (92)  is  found  to  be 

r  ca  2p  -[4-(l+v)2]n2  (93) 

4 

If  terms  in  n  and  higher  powers  of  n  are  neglected  in  the  definition 
of  r,  the  following  approximate  expression  results : 
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r  rj  1  - 


(9M 


Equations  (93)  and  (9]0  yield  the  approximate  buckling  stress  ratio 

p  «i{l+[2-(l+v)2]n2}  (99) 

which  is  vary  nearly  one-half  for  thin  shells. 

FF5  AND  FFh  CASES 

These  cases  of  edge  restraint  involve  the  replacement  of  the  condition 
Txv  =  ®  with  v  =  0  in  the  two  preceding  analyses;  the  second  row  of 
tnel  coefficient  matrix  in  equations  (64)  is  replaced  by 

( [l-Q(p2+vn2 )  ]  [l-Q(p2+vn2)  ]  [l-Q(p2+vn2 )  ]  [l-Q(p2+vn2)  ] }  (96) 

The  systematic  simplification  procedure  employed  in  the  two  previous  cases 
can  be  utilized  again.  The  algebra  is  slightly  more  involved,  but  the 
buckling  stress  ratios  for  both  FF3  and  FF4  restraints  are  one -half  plus 
small  correction  terms,  the  first  of  which  is  proportional  to  n^  as  in 
the  FF2  case. 

FC1  CASE 


The  simultaneous  satisfaction  of  equations  (56)  and  (60)  for  nontrivial 
values  of  A  through  is  assured  if  the  following  determinantal 

equation  is  satisfied: 


(97) 


If  the  following  operations  are  performed  on  the  above  determinant, 


1.  Multiply  the  columns  of  the  determinant  by  p  ,  p^,  p  ,  and  p  , 
respectively,  and  disregard  the  n2  factor  appearing  in  the  last 
three  rows. 

2.  Subtract  column  1  from  column  2,  and  subtract  column  3  from 
column  A. 

3.  Remove  the  factor  (p  -p,)(p  -p.)  =  1  from  the  determinant  as  in 
equation  (65). 
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4.  Rearrange  rows  and  columns  so  as  to  have  the  four  vanishing 
elements  occupy  the  upper  right  quadrant  of  the  determinant. 

Then  the  buckling  condition  becomes 


Q 

Q 

0 

0 

1 

1 

0 

0 

=  (Q-Q)[Q(p1+p2)-Q(p1+p2)]  =  0 

p2Q 

Tjl 

-Q 

-Q 

c. 

*1 

% 

(?1^2) 

(p-L+pJ 

(98) 

The  first  factor  in  equation  (98)  vanishes  for  p  =  1.  From  equation 
(78),  the  second  factor  vanishes  if 


p(r-Mn2+p)‘^2-(l-p2)'^2(r+4n2-p)‘^2  =  0  (99^ 

If  equation  (99)  is  multiplied  by  (r+4n2-p)^2  and  the  factor  (l-p2)1^2 
is  disregarded  as  in  equation  (8l),  the  following  result  is  obtained: 

r  =  2p  -  4n2  (100) 

If  both  sides  of  this  expression  are  squared  and  the  definition  of  r  is 
introduced,  the  quadratic  equation 

1  -  8r2p  +  16 =  (2p  -  4n2)2  (1C1N 

is  obtained;  this  result  may  also  be  written  as 

4p2  -  8n2p  -1=0  (102) 

The  positive  root  of  (102)  is 

p  =  ^  }j  1  +  4n^1  +  n2  (103) 

which  is  very  nearly  one-half  for  thin  shellE. 

FCg  CASE 


Equations  (56)  and  (6l)  will  be  satisfied  at 
determinant  in  equation  (97)  Is  replaced  by 


[P2r‘ 


p2Q 


PlQ 


] 


j=0  if  the  first  row  of 


(104) 
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If  a  sequence  of  elementary  operations  is  applied  to  the  determinant,  the 
buckling  condition  car.  be  simplified  to 


Q 

1 


P^Q 


Q 

1 


-2 


0 

0 

-Q 


o 

0 

-Q 


=  r(S-Q)2  =  0 


(105) 


The  factor  r  cannot  vanish  for  values  of  p  of  interest,  and  (Q-Q) 
vanishes  for  p  =  1;  therefore,  buckling  at  values  of  stress  less  than  the 
classical  value  is  not  possible.  A  discussion  of  buckling  at  p  =  1  is 
given  in  Appendix  III. 

FC3  CASE 


Equations  (56)  and  (62)  will  be  satisfied  at  |  =  0  for  nontrivial  values 
of  A1  through  A^  if 

Q  Q  Q  Q  I 


(1-QP^) 


(l-QPg)  (1-Qp^) 


(1-QP^) 


=  0 


P1 


P2 


(106) 


The  application  of  a  sequence  of  elementary  operations  reduces  the  buckling 
condition  to 

-(Q-Q)2  +  r.2[Q(p1+p2)-4(p1+p2)][Q(p1-+p2)-Q(p1+p2)]  ==  0  (107) 

If  equations  (7&),  (79),  and  (90)  are  utilized,  the  buckling  condition 
simplifies  to 


2(l-p2)  -  n2[r(2p2-l)  -  4n2  +  p]  =  0  (108) 

By  neglecting  terras  involving  n^  in  equation  (108)  and  using  equation 
(9U),  an  approximate  equation  is  obtained;  this  equation  is 

2(p+l)[l-p  +  n2(|  -  p)]  -  0  (109) 
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The  bracketed  term  in  this  expression  vanishes  (to  the  same  approximation 
in  terms  up  to  n^)  when 

2 

P  *  1  -  §-  (no) 

which  implies  buckling  at  slightly  les6  than  the  classical  value  of  stress . 


FC4  CASE 


The  FC4  case  represents  the  mosi  rigid  condition  of  tangential  restraint 
Since  buckling  stress  ratios  of  unity  and  slightly  less  than  unity  have 
been  found  for  the  FC2  and  FC3  conditions  of  restraint,  there  is  little 
doubt  that  buckling  will  occur  at  p  =  1  in  this  case.  It  is  indicated 
in  Appendix  III  that  p  =  1  is  always  a  possible  buckling  stress  ratio 
for  all  cases  of  edge  restraint;  this  fact,  along  with  the  intuitive 
physical  consequences  of  edge  fixity,  leads  to  the  conclusion  that  p 
should  be  unity  in  this  instance. 


SIMPLY  SUPPORTED  CASES 

The  SSI,  SS2,  and  SS3  cases  are  treated  in  detail  by  Hoff  and  the  author 
in  Reference  4.  The  buckling  stress  ratios  in  the  first  two  cases  are  the 
same  as  in  the  FC1  case,  while  the  SS3-supported  shell  can  buckle  only  at 
p  =  1  (this  i6  the  classical,  simply  supported  shell). 

The  same  argument  presented  for  the  FC4  case  applies  equally  well  to  the 
SS4  case.  Since  p  =  1  is  a  possible  buckling  stress  ratio  for  all  cases 
and  it  is  the  only  possible  value  for  buckling  in  the  SS3  case,  where  the 
tangential  restraint  is  less  severe,  one  must  conclude  that  p  =  1  i6  the 
buckling  condition  for  SS4-aupported  shells. 

COMPLETELY  CLAMPED  CASES 

All  the  cases,  CC1  through  CC4,  were  solved  numerically  by  Hoff  and 
Soong  [6].  An  examination  of  their  results  irdicates  that  if 

5>10S  :m) 

where  L  is  the  shell  length,  then  p  is  essentially  unity  for  all  clamped 
edge  conditions. 

The  results  of  Reference  6  can  be  verified  readily,  and  the  next  section  is 
devoted  to  a  short  analysis  of  the  CC1  case. 
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CC1  CA.jE 


The  nontrivial  satisfaction  of  equations  (ljtt  and  (60)  at  l  =  0  is 
assured  if  the  following  determinantal  equation  is  satisfied: 


1 

Q 

Pi 

Pi  Q 


l 

Q 

P2 

PgQ 


1 

Q 

Pi 

PiQ 


1 

Q 


P? 

p2Q 


=  o 


(112) 


The  rows  in  the  above  determinant  reflect  the  vanishing  of  w,  <r 

t°  ,  respectively.  X1 

xy>1 


and 


,  , 


If  column  1  is  subtracted  from  column_2,  column  3  is  subtracted  from 
column  U,  the  factor  (pg  -  P]_)(pg  -  p]_)  =1  is  noticed,  and  the  columns 
are  rearranged  in  the  above  determinant,  equation  (112)  can  be  replaced  by 


110 
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Q  Q  0 

Pl  ?!  1 

PjQ  plQ  Q 


o 

l 

Q 


■(Q-Q)2  =  0 


(113) 


(Q-Q)  vanishes  only  at  o  =  1  ,  so  buckling  cannot  occur  at  values  of 

stress  less  than  the  classical  value. 


CC2,  CC3,  AND  CCU  CASES 

Each  of  these  cases  implies  more  rigid  edge  fixation  than  the  CCl  case; 
it  is  realized,  therefore,  that  p  =  1  corresponds  to  buckling  in  these 
cases. 
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CONCLUDING  REMARKS 


A  theoretical  analysis,  based  upon  linear  theory,  of  the  buckling  behavior 
of  axially  compressed  circular  cylindrical  shells  has  been  undertaken  to 
determine  the  effect  of  edge  restraint  upon  the  buckling  stress;  and  over¬ 
all  evaluation  of  this  effect  has  been  made  by  studying  the  l6  limiting 
cases  of  restraint  which  correspond  to  the  vanishing  of  generalized 
forces  and  generalized  displacements  at  an  edge.  Of  the  l6  limiting 
cases,  8  permit  buckling  to  occur  at  stresses  less  than  the  classical 
value  of  buckling  stress. 


The  buckling  stress  ratios  for  all  the  cases  are  displayed  in  the  table  on 
page  38.  While  buckling  in  one  case  occurs  at  8  value  of  stress  only 
slightly  less  than  the  classical  one,  buckling  stresses  approaching  one- 
half  of  this  value  are  predicted  in  six  instances  which,  interestingly 
enough,  represent  a  rather  wide  variety  of  restraining  conditions;  a 
further  1?  percent  stress  reduction  has  Deen  found  previously  by  Nachbar 
and  Hoff  for  completely  unrestrained  cylinders  of  sufficient  length.  In 
addition  to  the  seven  cases  of  restraint  which  have  been  analyzed  hero  for 
the  first  time,  the  work  of  previous  investigators  for  the  remaining  cases 
has  been  reestablished  and  supplemented  as  well.  Since  the  effect  of  shell 
length  has  been  shown  to  be  of  minor  importance  for  shells  of  practical 
interest,  the  analysis  is  simplified  by  treating  Bern! -inf inite  shells,  and 
this  simplification  permits  solutions  to  be  obtained  in  essentially  closed 
form. 

On  the  basis  of  the  results  obtained,  both  by  previous  investigators  and 
by  the  present  analysis,  it  is  concluded  that  the  mode  of  edge  restraint 
can  have  a  very  significant  effect  upon  the  buckling  stress  of  thin, 
axially  compressed  circular  cylindrical  shells.  This  effect  has  been 
widely  overlooked  in  the  past,  a  possible  reason  being  the  absence  of  a 
pronounced  shell  length  effect  upon  the  buckling  stress  in  both  theory  and 
experiment  except  when  extremely  short  shells  are  considered. 
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APPENDIX  I 


PROOF  M  THE  ROOTS  p  THROUGH  p^  HAVE  NEGATIVE  REAL  PAR^S  WHEN  p  <  1 

The  roots  p1  through  will  have  negative  real  parts  if  the  T>a Trimeters 
and  a2  xdefined  in  equr+ion  (44)  are  negative;  these  parameters  are 


ax  =  -(1-p)1/2  -(r+4n2-p)l/2 

(114) 

a2  =  (l-p)1//2  -  (r+4n2-p)1/2 

(115' 

If  p  <  1,  Ct^  will  be  negative  if 

r  +  4r2  -  p  >  0  (116) 

By  the  definition  inequation  (4l)  for  r,  equation  (116)  can  be  written 

[(4n2-p)2  +  l-p2]1/2  +  (4n2-p)  >  0,  (117) 

which  proves  that  is  obviously  negative  for  p  <  1. 

Equation  (115)  implies  that  GC^  will  be  negative  If 

(r+4n2-p)1/2  >  (1-p)1/2  (116) 

This  is  equivalent  to  the  requirement  that 

r  >  1  -  n2  (119) 

Now  r  can  be  conveniently  written  as 

r  =  [(4n2-p)2  +  1-p2]1/2  =  [  (l-4n2)2  ♦  6n2(l-p)]l/2  (120) 

Thus,  if  equation  (L20  )  is  substituted  into  equation  (119),  It  is  clear 
that  a2  is  also  negative  for  p  <  1. 

These  conclusions  are  stated  in  Reference  14,  but  a  proof  is  not  included. 
/  A  less  direct  proof  is  also  given  in  Reference  4. 
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APPENDIX  II 


\ 


IDENTITIES  INVOLVING  THE  ROOTS  pg,  J>y  AMD  p^ 

Equations  (28)  through  (4l)  define  and  describe  the  roots  of  equation  (24) , 
The  four  roots  of  interest,  which  possess  negative  real  parts  for  p  <  1, 
can  be  written  as 


Pl  =  |[-Q1/2  -  (Q^n2)1/2] 

P2  =  ICq1/2  ”  (Q+4n2>l/2] 


?1  =  |t-Ql/2  -  (5+4n2)l/2]  =  p. 


p0  =  i[*l/2  -  (Q^n2)1/2)  - 


P2  =  2 

These  roots  satisfy  the  following  equations: 

.2  2 


p4 


(121) 

(122) 

(123) 

(124) 


2  2 

P1  0l/2 

Pl  ' 


P2-n 


,1/2 


-2  2 

Pl-n 


Ql/2, 


-2  2 

P2  _  nl/2 


(125-128) 


On  the  basis  of  equations  (121)  through  (124),  it  is  clear  that 

p2  ’  P1  =  ^ 2  P1  +  p2  =  -(Q^2)1^2 
p2  -  Pi  -  5l/a  ?!  *  P2  -  -<5^na)1/2 

(p2*P1)(p2-P1)  »  -  1 

Equations  (125)  through  (126)  provide  the  identities 


(129-133) 


+  P2  =  2n2  +  (p2-P;L)Ql/2  =  2n2  +  Q 

(134) 

+  p|  =  2n2  +  (Pg-P!^2  =  2n2  +  Q 

(135) 
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With  the  aid  of  equations  (40)  and  (4l),  It  is  clear  that 

(Px  +  P2)(Pi  +P2^  =  l^+4n2)i'/2!  =  r  ( 1 36) 


The  Identities  given  here,  in  addition  to  those  derived  in  the  development 
of  the  solutions,  are  sufficient  for  the  derivation  of  the  buckling 
solutions.  Some  of  the  above  Identities  ware  originally  given  by  Nachbar 
in  Reference  x4,  and  others  can  be  found  in  References  4  and  6. 
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APPENDIX  III 


PR  PC?  THAI  p  =  1  IS  ALWAYS  A  BUCKLING  STRESS  RATIO  FOR  LONG  SHELLS 


By  assigning  the  value  unity  to  the  stress  ratio  p,  equations  (28),  (29), 
(3l),  and  (3u)  yield  the  following  simple  expressions: 

Q  -  Q  =  -1  (137) 

«1/2  =  1  (138) 

Q1/2^-i  (139) 


If  these  results  are  substituted  into  equations  (l2l)  through  (124),  it  is 
found  that 


Px  =  -  |[1  +  (l-^n2)1/2] 

P2  -  |[1  -(l-4r2)1//2] 

P3  =  P2 

p4  =  P1 


fi4o) 

<  l4l) 

(142) 

(143) 


Equations  (137),  (192),  and  (143)  imply  that  the  first  and  fourth  columns, 
as  well  as  the  second  and  third  columns,  of  all  the  buckling  determinants 
written  for  semi-infinite  shells  become  identical  when  p  =  1.  Therefore, 
the  vanishing  of  all  these  determinants  at  p  =  1  is  assured. 


The  question  of  whether  buckling  modes  for  p  =  1  can  be  found  to  satisfy 
a  given  set  of  edge  restraint  conditions  is  rather  academic.  It  is 
sufficient  merely  to  refer  to  the  numerical  results  obtained  by  Hoff  and 
Soong  in  Reference  6  for  shells  of  finite  length.  These  results  indicate 
that  even  for  a  CC4 -supported  shell  of  finite  length,  p  is  essentially 
unity  whenever 


L 

R 


t 

>  10  R 


where  L  is  the  length  of  the  shell.  If  t/R  =  l/lOO,  then  a  shell  must 
be  shorter  than  its  radius  before  the  buckling  stress  is  noticeably 
increased  above  the  classical  value.  Since  the  CC4  case  represents  the 
most  severe  condition  of  edge  restraint,  unity  is  an  upper  bound  for  uhe 
buckling  stress  ratio  which  can  be  applied  to  all  shells  whose  lengths 
exceed  their  respective  radii,  regardless  of  the  condition  of  edge 
restraint . 
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A  broad  analytical  study  of  the  effect  of  edge  restraint  on  the  buckling  stress  of 
geometrically  perfect,  axially  compressed  circular  cylindrical  shells,  based  on 
linear  elastic  theory,  is  presented.  An  overall  evaluation  of  this  effect  is  made  by 
studying  the  16  limiting  cases  of  restraint  which  correspond  to  the  vanishing  of 
generalized  forces  and  generalized  displacements  at  an  edge.  Of  the  16  limiting 
cases,  7  of  which  are  studied  for  the  first  time,  it  is  found  that  8  permit  buckling 
to  occur  at  stresses  less  than  the  classical  theoretical  value  of  buckling  stress. 

The  fact  that  both  theory  and  experiments  indicate  that  the  effect  of  shell  length 
on  the  buckling  stress  can  usually  be  disregarded  is  used  in  the  analysis.  This 
fact  suggests  that  shells  may  be  assumed  to  be  semi-infinite;  as  a  consequence, 
formulas  for  buckling  stresses  are  determined  in  essentially  closed  form. 


